Multi-color Shallow Decay and Chromatic Breaks in the 
GRB 050319 Optical Afterglow 



K.Y. Huang 1 , Y. Urata 2 - 3 , P.H. Kuo 1 , W.H. Ip 1 , K. Ioka 4 , T. Aoki 5 , C.W. Chen 1 , 
W.P. Chen 1 , M. Isogai 5 , H.C. Lin 1 , K. Makishima 3 ' 6 , H. Mito 5 , T. Miyata 5 , Y. 
Nakada 5 , S. Nishiura 7 , K. Onda 2 , Y. Qiu 8 , T. Soyano 5 , T. Tamagawa 3 , K. 
Tarusawa 5 , M. Tashiro 2 , and T. Yoshioka 9 

ABSTRACT 

Mult i- wavelength B, V, R, I observations of the optical afterglow of 
GRB 050319 were performed by the 1.05-m telescope at Kiso Observatory and 
the 1.0-m telescope at Lulin Observatory from 1.31 hours to 9.92 hours after 
the burst. Our R band lightcurves, combined with other published data, can be 
described by the smooth broken power-law function, with a± = —0.84 ±0.02 to 
a2 = — 0.48±0.03, 0.04 days after the GRB. The optical lightcurves are char- 
acterized by shallow decays — as was also observed in the X-rays — which may 
have a similar origin, related to energy injection. However, our observations in- 
dicate that there is still a puzzle concerning the chromatic breaks in the R band 
lightcurve (at 0.04 days) and the X-ray lightcurve (at 0.004 days) that remains 
to be solved. 
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1. Introduction 

The GRB afterglow as perceived in the X-ray, optical and radio wavelengths is now 
understood to be the result of the collision between relativistic ejecta from the gamma-ray 
bursts and the interstellar medium (ISM). A comparison of afterglow lightcurves obtained at 
different wavelengths gives important information about the surrounding ISM environment 
and the interaction processes. Such analyses can also provide essential input for theoretical 
models. Recently, the pace of this type of activity has quickened significantly, stimulated 
by the capabilities of the quick response and accurate localization of the GRB by the Swift 
Satellite (Gehrels et al. 2004). This has meant that the number of GRB optical afterglow 
detections in the first several hours by ground-based telescopes has recently increased signif- 
icantly. It is interesting to note that the observations by Swift of the early X-ray emissions 
from a number of GRBs reveal a canonical behavior. The X-ray lightcurves can be divided 
into three distinct power law segments (Nousek et al. 2006). Some X-ray and optical ob- 
servations show that the evolution of both lightcurves changes at the same time (Blustin 
et al. 2006; Rykoff et al. 2006), however chromatic breaks were also found in some cases 
(Fan & Piran 2006; Panaitescu et al. 2006). The nature of the afterglow early breaks in 
the lightcurves is thus uncertain. A detailed comparison of changes in the evolution of the 
optical, radio and X-ray lightcurves should therefore be very interesting. This kind of physi- 
cal study demands both a well-coordinated observational program and careful data analysis. 
We use GRB 050319 which has comprehensive observational coverage in both the X-ray and 
optical wavelengths may be used as just such an example. 

GRB 050319 was detected by the Burst Alert Telescope (BAT) instrument onboard 
the Swift satellite on March 19, 2005 at 09:31:18.44 UT (Krimm et al. 2005). However, 
a re-analysis of the BAT data showed two flares, which indicated that GRB 050319 had 
already started 137 sec before the trigger. The 15-350 keV fluence for the entire burst 
duration of Tgo = 149.6 ± 0.7 sec has been estimated to be 1.6 xlO -6 erg cm -2 (Cusumano 
et al. 2006). The X-ray emission of GRB 050319 after the burst was monitored by the 
XRT from 225 sec to 28 dayfl Two breaks in the emission curves were found (Cusumano 
et al. 2006). The initial sharp decline can be described by a power-law with an index of 
cti = —5.53 ± 0.67 to be followed by a<i = —0.54 ± 0.04 after 0.004 days since the burst. 
The unusually flat decline in the second part might have been caused by continuous energy 
injection (Cusumano et al. 2006). At about 0.313 days after the burst, the power-law index 
changed to a 3 = —1.14 ± 0.2 which can be readily explained as a jet break or a reduction in 
the energy injection (Cusumano et al. 2006; Zhang et al. 2006) . 



1 The burst time in the article is 09:29:01.44 UT, 137 s before the BAT trigger 
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The early optical afterglow emission at 230 s was observed by the UVOT telescope on 
Swift (Mason et al. 2006) and by two ground-based robotic telescopes, ROTSE-III (Quimby 
et al., 2006) and RAPTOR (Wozniak et al. 2005). The best single power-law fit of unfiltered 
data from ROTSE-III and RAPTOR indicates that a = -0.854 ± 0.014. A number of 
optical observatories have joined the follow-up observations (see Yoshioka et al. 2005; Torii 
et al. 2005; Sharapov et al. 2005a,b; George et al. 2006; Misra et al. 2005; Kiziloglu et al. 
2005; Greco et al. 2005). The spectral measurements of the afterglow by the Nordic Optical 
Telescope (NOT) indicate that was redshift z = 3.24 of this event (Fynbo et al. 2005). 



2. Observations and Analysis 

After receiving the GRB alert message from Swift and the afterglow position was re- 
ported by Rykoff et al. (2005), the Target-of-Opportunity procedures of the East- Asia GRB 
Follow-up Observation Network (EAFO^@, Urata et al. 2005) were immediately carried out. 
A series of multi-band follow-up observations were successfully performed by the 1.0 5-m 
Schmidt telescope of the Kiso Observatory in Japan and the Lulin One-meter Telescope 
(LOT) in Taiwan. Photometric B and R images were obtained at the Kiso site with a 2k 
x 2k CCD camera (Urata et al. 2005) between 0.055 and 0.326 days after the burst. A 
number of parallel B, V, R and I images were obtained by LOT with a PI1300 CCD camera 
(Kinoshita et al. 2005) from 0.080 to 0.413 days after the burst. 

The standard routine including bias subtraction and dark subtraction; flat-field correc- 
tions were employed with the appropriate calibration data needed to process the data using 
IRAF. The afterglow can be clearly seen in the images. The signal-to-noise ratio was im- 
proved by combing the LOT B band data with median filtering. The DAOPHOT package 
(Stetson 1989) was then used to perform point-spread function (PSF) fitting for the GRB 
images. Four field stars were used to create a PSF model which was applied to the optical 
afterglow of each GRB image. For absolute photometric calibration, we used calibrated data 
of the GRB field obtained by the USNOFS 1.0-m telescope (Henden 2005). The photometric 
error and the systematic calibration error were included in the magnitude error estimation^]. 



2 EAFON web-page: |http: / /cosmic.r iken.jp/grb/eafon/ 

3 The errors in this article were quoted for 68% (1-cr) confidence level 
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3. 



Results 



3.1. Lightcurve 



Figure 1 shows the multi-band lightcurves of the GRB 050319 afterglow. Besides our B, 
V, R and I band data (Table 1) we also included the R band measurements from ROTSE-III 
(Quimby et al. 2006), RAPTOR (Wozniak et al. 2005) and several GCN reports (Greco et 
al. 2005; Kiziloglu et al. 2005; Misra et al. 2005; Sharapov et al. 2005a,b). In addition, 
we also made use of several B and V band measurements taken by the Swift UVOT (Mason 
et al. 2006). The GCN R band points were re-calibrated using the GRB 050319 field stars 
reported by Henden (2005) so they could be plotted on the same magnitude scale. The 
magnitude differences between photometric field stars in Henden (2005) and in USNO-A2.0, 
USNO-B1.0 stars are +0.18 and —0.22 mag, respectively. We re-measured the reference stars 
from Greco et al. (2005) from the LOT R band images and obtained the average magnitudes 
and rms errors. 

After fitting the B,V,R and I band lightcurves to a single power-law expression F oc t a , 
where a is the index and t is the time after the burst, we get a = —0.56 ± 0.06 {x 2 /v= 2.90 
for v = 19) for the B band, a = -0.65 ± 0.03 {ydjv= 2.60 for v = 27) for the V band, 
a = -0.59±0.01 (x 2 /z/= 5.3 for v = 97) for the R band, and a = -0.52±0.15 (x 2 /v= 7.7 for 
v — 9) for the I band. This single power-law fitting indicates that these lightcurves, obtained 
with different filters have similar power-law decay even though the reduced chi-square values 
are relatively large. 

Since the data sets of the V and R measurements are more complete, it is possible with 
following expression to attempt the fitting of the corresponding lightcurves with a smoothly 
broken power-law function: 



where % is the break time, ot\ and a>2 are the power-law indices before and after F u ^ is 
flux at break t b , and k is a smoothness factor. For the V band, we obtain a\ = —0.87 ±0.21, 
a 2 = -0.49 ± 0.05, t h = 0.042 ± 0.058 days and k = -30 (x 2 /^= lA % v = 24). For 
the R band, we obtain a x = -0.84 ± 0.02, a 2 = -0.48 ± 0.03, t h = 0.046 ± 0.008 days, and 
k = —21 (x 2 /^= 2.24 for v = 90). This result implies a mild break in both the V and the 
R band lightcurves at around 0.04 days after the occurrence of the GRB. 

Taking t b = 0.04 days, we fit the data in the B and I bands to a respective power-law 
before and after the break. In this manner, we find a,\ = —0.79 ± 0.09 {y 2 jv = 1.09 for 
v = 7); a 2 = -0.36 ± 0.05 {yfjv = 1.23 for v = 9) for the B band and a 2 = -0.52 ± 0.15 



F(u,t) 



[(t/tb)- fcai + (t/*b)- fca2 ] (1/fc) ' 



2 1/k F u . 



(1) 
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(x 2 /V = 7.7 for v — 9) for the I band. The best-fit parameters for the B, V, R, and / bands 
are summarized in Table 2. Our results show not only the clear presence of mild breaks in 
the V, R band lightcurves but a flattening trend after the break. Furthermore, our R band 
slope before the break (a ~ —0.84) is in agreement with the corresponding value derived by 
Quimby et al. (2006) for the interval between 0.0019 and 0.05 days after the burst. 

3.2. Color and spectral flux distribution 

Our multi-wavelength observations indicate that median colors between 0.07 and 0.35 
days are V-R = 0.45±0.11, R-I = 0.46±0.10, and B-V = 0.84±0.14. These values have 
been corrected for foreground reddening of E(B — V) = 0.011 mag (Schlegel et al. 1998). 
The V — R and R — I colors so derived are consistent with those of the typical long GRBs 
(Simon et al. 2001), but the B — V color is slightly redder than those of the typical long 
GRBs (B — V — 0.47 ± 0.17). The larger B — V value may imply a certain absorption effect 
because the redshift of GRB 050319 was determined to be 3.24 (Fynbo et al. 2005). 

The B, V, R and I magnitudes have been further converted to fluxes using the effec- 
tive wavelengths and normalizations of Fukugita et al. (1995). The effect of the Galactic 
interstellar extinction has been corrected. Figure 2 shows two samples of spectral energy 
distribution obtained by LOT at 0.13 and 0.21 days after the occurrence of GRB 050319. A 
drop in the B band flux at about 4380A can be clearly seen. We subsequently fitted the flux 
distribution of V, R and / bands with a power-law function F(u, t) oc u 13 ; here F(u, t) is the 
flux at frequency v with a certain t and (3 is the spectral index. We find that (3 = — 1.08±0.05 
( x 2 /is= 0.05 for v = 1) at 0.13 days and that (3 = -1.08 ± 0.32 {^jv = 2.3 for u = 1) at 
0.21 days. Our result (f3 = —1.08 with a rms error 0.23) is consistent with the X-ray fitting 
value ((3 = -0.69 ± 0.06) in a 3-a level. 

With a redshift of 3.24, the Lya absorption feature would shift into the B bandpass, 
causing reduction of the afterglow flux in the B band. To correct for this absorption effect, 
we used the formulation derived by Yoshii et al. (1994), in which the optical depth is a 
function of the observed wavelength and source redshift. With the computed optical depth 
in the B band, and a spectral slope of (3 — —1.08, we found the expected B band magnitude 
after Lya absorption to be 21.33 ±0.05 at 0.13 days. This value compares very well with our 
observed value of B = 21.26 ± 0.17 at 0.13 days after correction for Galactic extinction. The 
drop at the B band is hence fully produced by the Lja absorption and no spectral breaks 
should have taken place during our observation. 
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4. Discussion and Summary 

It is important to note that Swift found two breaks at 0.004 and 0.313 days after the 
burst in the X-ray afterglow observations (Cusumano et al. 2006), but we only found a single 
break in our V and R lightcurves (see Figure 1). In the following since there are more data 
points available for the R band data we will focus on this. It is useful to remember that a.\ 
= —0.84 and a 2 = —0.48 at the break time of t h = 0.04 days. 

4.1. Before the optical break ( t < 0.04 days) 

The slope ol\ (=—0.84) is consistent with the typical range of a — —0.62 to —2.3 for 
many well observed GRBs. According to the standard afterglow model relating the power- 
law index (a) to the power-law index (p) of an electron spectrum (Zhang & Meszaros, 2004; 
Dai & Cheng, 2001), the corresponding value for a\ = —0.84 is p — 2.1, which is in agreement 
with the constant-density ISM model with slow cooling in which p > 2 for v m < u opt < u c 
(v m is the typical frequency; u opt is the optical frequency; u c is the cooling frequency). In 
light of the XRT observations, the first break was likely caused by the transition from the 
tail end of the low energy prompt emission to the afterglow phase (Zhang et al. 2006). 
However, it is important to note that the X-ray break at 0.004 days (where the steeper slope 
becomes shallow) is not accompanied by an R band break. At the same time, the power-law 
decay slope in the X-ray (~ — 5.53) and the R band (~ — 0.84) are quite different. This 
is an indication not only that the behaviors of the X-ray afterglow and optical afterglow 
of the GRB 050319 event are different, but also suggested that the afterglow phase already 
dominated the optical bands when the optical emission was first detected. 

4.2. Shallow Decay 

The power law index becomes shallow after the break (t b = 0.04 days). Neither the jet 
(Rhoads 1999) nor the break frequencies across the optical wavelength (Sari et al., 1998) 
suitably explain the break they see in the GRB 050319 lightcurves. As discussed before, 
the X-ray lightcurve between the two breaks 0.004 and 0.313 days after the burst is also 
characterized by shallow decay. Zhang et al. (2006) suggested that such behavior is related 
to continuous energy injection into the ISM. For a long-lasting central engine, the energy 
injection rate is E(t) oc t~ q and with q < 1 (Zhang & Meszaros 2001). For slow cooling 
in the ISM, the temporal index can be expressed as : a = [(2p — 6) + (p + 3)g]/4 = 
[(q — 1) + (2 + q)/3}/2, when u m < v < u c . Using this formulation, Zhang et al. (2006) 
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obtained q = 0.6 and p = 2.4 from the X-ray observations. With a = —0.48 and (3 = —1.08 
from the R band observations, we find that q = 0.72 and p = 2.12. The results not only 
indicate that the electron spectrum power law index is the same before and after the break, 
they also compare well with the results of Zhang et al. (2006). These results indicate that 
the shallow decays evidenced by both X-ray and optical afterglows could be of a similar 
origin, related to a continuous energy injection mechanism. 

According to the energy injection model, would also expect an X-ray break at the time 
of the optical flattening break, because the onset of the energy injection should also alter 
the X-ray temporal index. However, such an X-ray break is not observed, which suggests 
that some modifications to the injection model may be needed. As mentioned in § 4.1, the 
X-ray break at 0.004 days was not accompanied by a break in the R lightcurve. Although 
a chromatic break in the X-ray was found at 0.004 days and in the optical region at 0.04 
days, the lightcurves at both wavelengths indeed showed shallow decay after the breaks, 
which can be explained by the energy injection model. However, it is difficult for energy 
injection from 0.004 days to 0.04 days to affect only high energies. This difficulty indicates 
that energy injection is an imperfect mechanism for explaining the shallow optical or X-ray 
phase associated with the GRB 050319 event. 

Several models have recently been proposed to explain the shallow decay effect. Using 
the multiple-subjet model (Nakamura, 2000), Toma et al. (2006) invoked the superposition 
of afterglows from many off- axis subjets. Eichler & Granot (2006) favored a combination 
of the tail of prompt emission model with the afterglow emissions observed from a viewing 
angle outside the edge of the jet. These arguments hence suggest that the multiple-subjet 
model and the patchy-shell model (Kumar & Piran, 2000) might provide a theoretical basis 
for explaining the observed shallow decays in the X-ray and optical lightcurves. It is inter- 
esting to note that in order to sustain the shallow decay process these models all require 
high gamma-ray efficiency (75-90%); additional mechanisms such as prior activity (Ioka et 
al. 2005) and time-dependent shock generation (Fan & Piran, 2006) have been also pro- 
posed. Comprehensive multi-wavelength observations, such as those reported here, provide 
important keys to improve these models. 

Finally, the second break in the X-ray emissions (~ 0.313 days after the burst), has 
been interpreted as being due to an unusual flat jet break (Cusumano et al. 2006), Zhang 
et al. (2006) however provided an alternate explanation, a sudden cessation of the energy 
injection. In both interpretations, the corresponding break should appear in both the X-ray 
and optical lightcurves. This effect cannot be clearly identified in our measurements until 
0.413 days after the burst. The lack of data for the subsequent time interval could lead to 
uncertainty in the power-law fitting. We thus cannot fully exclude the existence of a second 
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break in the optical lightcurves. However, Panaitescu et al. (2006) have studied several 
afterglows. They found the shallow power-law decay evidenced by the X-ray emissions to 
steepen about 0.04—0.17 days after the burst, although there was no accompanying break 
found in the optical range. They suggest that such chromatic X-ray breaks may be common. 
The chromatic breaks (e.g. the shallow X-ray phase becomes steeper, with no accompanying 
optical break) may be caused by differences in the X-ray and optical outflow (Panaitescu et al. 
2006), or by changes in the typical electron energy parameters (the so-called microphysical 
parameters) at the end of energy injection (Panaitescu 2006). 

In summary, our analysis of the optical multi- wavelength observations of GRB 050319 
made at Kiso and Lulin compared with the X-ray observations from Swift found the following 
major results: 

• The B, V, R, and / band lightcurves displayed unusual shallow decays. 

• The R lightcurve can be described by a smooth broken power-law function. a.\ ~ —0.84 
becomes shallow (a 2 ~ —0.48) 0.04 days after the occurrence of the GRB. 

• The shallow decay observed in the X-ray and optical lightcurves may have a similar 
origin related to energy injection. However, our observations indicate that there is still 
a major puzzle remaining concerning the chromatic breaks in the R band lightcurve 
(at 0.04 days) and the X-ray lightcurve (at 0.004 days). 

• Our calculations revealed that the drop in spectral energy distribution was fully caused 
by a shift in the Ly« absorption to the B bandpass at z = 3.24. 
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Table 1: GRB 050319 Optical Afterglow Photometry 



Days after GRB b 


Filter 


Magnitude 


Site 


0.05443 


R 


19.01±0.08 


Kiso 


0.06023 


B 


20.42±0.13 


Kiso 


0.08037 


V 


19.75±0.07 


Lulin 


0.10439 


I 


19.09±0.07 


Lulin 


0.27821 


B 


20.92±0.12 


Kiso 


0.33424 


I 


19.91±0.09 


Lulin 


0.40097 


V 


20.68±0.12 


Lulin 


0.41346 


R 


20.20±0.07 


Lulin 



a Table 1 is published in its entirety in the electronic edition of the Astrophysical Journal. 

b Thc burst time is 2005 March 19, UT 09:29:01.44 

c The magnitudes are not corrected for Galactic extinction. 
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Table 2: Fitting results of the GRB 050319 lightcurves : The B and / band data were fitted 
by a respective power-law model [F oc t a ] before and after the break (tj, = 0.04 days). On 
the other hand, the V and R data were fitted by a smoothly broken power-law [equation 

(1)1- 



Filter 




ot 2 


h (day) 


B 


-0.79±0.09 


-0.36±0.05 




V 


-0.87±0.21 


-0.49±0.05 


0.042±0.058 


R 


-0.84±0.02 


-0.48±0.03 


0.046±0.008 


I 




-0.52±0.15 
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Days after the burst 

Fig. 1. — Optical lightcurves of GRB 050319 : The solid lines present the best fit by the 
single power-law model [F oc t a ] for the B band (a = —0.56 ± 0.06) and the I band 
(a = —0.52 ± 0.15). The dashed line indicates the best fit by the smooth broken power-law 
[equation (1)] with the V band (on = -0.87 ± 0.21, a 2 = -0.49 ± 0.05, t h = 0.042 ± 0.058 
days) and the R band (a x = -0.84 ± 0.02, a 2 = -0.48 ± 0.03, t b = 0.046 ± 0.008 days). 
These observations are based on our EAFON data and Quimby et al. (2006), Wozniak et 
al. (2005), Mason et al. (2005), Greco et al. (2005), Kiziloglu et al. (2005), Misra et al. 
(2005), and Sharapov et al. (2005a,b). The dotted lines represent the break times of X-ray 
afterglows at 0.004 and 0.31 days after the burst. The two breaks were not found in X-ray 
observations, but a mild break seems to exist in the V and R band lightcurves. 
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Fig. 2. — Spectral energy distribution of GRB 050319 between the V, R and I bands at 0.13 
and 0.21 days after the burst (corrected for Galactic extinction). The solid and dashed lines 
indicate the best fit by the power- law model [F{y) oc v% where (3 = -1.08 ± 0.05 for 0.13 
days and j3 = —1.08 ± 0.32 for 0.21 days, respectively. 



